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SOME  DATA  ON  LARGE-SCALE  FIELD  CHARACTERISTICS  OF  HORIZONTAL  VELOCITY  COMPONENTS 

IN  THE  OCEAN 

[Ozmldov,  R.  V. , Nckotoryye  dannyye  o krupnomasshtabnykh  kharakteristikakh  polya 

gor izontal ' nykh  komponent  skorosti  v okeane,  Izvestiya  AN  SSSR.  Seriya  geof Izicheskaya, 

No.  11,  1964,  pp.  1708-1719;  Russian] 


^A  series  of  statistical  field  characteristics  of  horizontal  velocity  /170£ 

components  In  the  ocean  are  calculated  by  using  the  method  of  analytical 
filtering  of  the  investigated  functions  in  different  filter  passbands.  It 
It  shown  that  the  spectral  energy  density  function  of  velocity  fluctuations 
Is  In  good  agreement  with  the  5/3  power  law  up  to  24-h  fluctuation  periods. 

The  values  of  a series  of  terms  in  the  equation  of  motion  of  ocean  waters 
are  also  estimated.. 


At  the  present  time,  onAof  the  important  problems  of  naval  hydrodynamics  is 
the  determination  of  the  role\pf  turbulent  friction  in  liquid  motion  in  the  coastal 
and  open  parts  of  the  ocean,  jhis  problem  may  be  solved  in  two  different  ways.  The 

first  consists  in  the  direct  calculation  of  the  derivatives  of  turbulent  stresses  - 

Pujuj  in  the  Reynolds  equations,  which  describe  the  turbulent  flow  studied: 


dP  . d(—  puj'u/y 

dxi  dxj 


(1) 


where  1,  J ■ 1,  2,  3 (the  repeating  index  denotes  summation);  x^  are  Cartesian  coor- 
dinates; Uj  are  the  averaged  velocity  components;  are  the  fluctuational  or  tur- 
bulent velocity  components,  equal  to  the  deviations  of  the  corresponding  instantaneous 
(nonaveraged)  velocity  components  from  their  averaged  values;  F^  are  the  components 
of  external  forces;  p is  the  pressure;  p is  the  density  of  the  liquid;  and  t is  the 
time.  The  upper  bar  in  the  Reynolds  equations  denotes  a given  type  of  averaging 
(for  example,  time  averaging)  with  a definite  scale  (period).  Terms  describing  the 
action  of  molecular  viscosity  forces  were  omitted  from  Eqs.  (1),  since  they  are 
negligibly  small. 

The  second  method  of  studying  the  role  of  turbulent  friction  in  the  sea  involves 
the  introduction  of  certain  hypotheses  relating  turbulent  stresses  in  the  Reynolds 
equations  to  the  field  of  averaged  velocities  in  the  flow.  The  various  coefficients 
thus  introduced  cannot  be  determined  theoretically,  but  can  be  estimated  experimen- 
tally in  the  turbulent  current  studied.  However,  in  most  cases  these  coefficients 
are  not  universal  constants,  but  are  themselves  substantially  dependent  on  the  char- 
acteristics of  the  flow  velocity  field  and  many  other  factors.  In  addition,  the 
hypotheses  themselves  are  frequently  inadequately  substantiated  and  their  applica- 
bility may  be  called  into  question.  Therefore,  the  first  method  of  studying  tur- 
bulent friction  in  the  sea,  not  Involving  the  introduction  of  any  hypotheses  or 
coefficients,  but  based  on  a direct  calculation  of  the  corresponding  terms  in  the  /170! 
Reynolds  equations,  appears  to  be  preferable. 


Numbers  in  the  right  margin  indicate  pagination  in  the  orLginal  text. 


A major  objective  In  this  case  Is  not  only  to  determine  the  values  of  turbulent 
stresses  and  their  derivatives,  but,  most  importantly,  to  study  their  dependence  on 
the  averaging  scale  Tq,  used  in  resolving  the  velocity  field  Into  the  averaged 
(regular)  and  f luctuatlonal  (turbulent)  components.  The  special  Importance  of  this 
question  for  sea  conditions  Is  due  to  the  presence  in  sea  turbulence  of  eddies  of 
the  most  diverse  dimensions;  this  must  inevitably  result  in  a substantial  dependence 
of  the  averaged  quantities  and  fluctuation  velocities  (and  hence,  Reynolds  stresses) 
on  the  scale  Tq  used  in  a given  problem. 

To  calculate  the  forces  of  turbulent  friction  from  experimental  data,  It  Is 
necessary  to  have  a long  series  of  observations  of  current  velocities  at  a series 
of  ocean  points  sufficiently  close  to  each  other.  Unfortunately,  no  such  data  are 
available  to  the  investigators  at  the  present  time.  However,  a sufficient  number 
of  day-to-day  recordings  of  horizontal  velocity  components  [illegible]  = u and 
u^  = v at  individual  points  of  the  ocean  are  now  available,  which  make  it  possible 

to  calculate  the  Reynolds  stresses  -pu'v'  = -pv'u',  -pu'2,  -pv'^  and  to  study  the 
dependence  of  these  stresses  on  the  averaging  period  Tq. 

At  the  present  time,  most  recordings  of  current  velocities  in  the  sea  are  ob- 
tained by  means  of  BPV-2  type  printing  meters.  By  virtue  of  their  structural  charac- 
teristics, BPV-2  current  meters  placed  in  the  sea  record  instantaneous  directions  of 
the  horizontal  current  velocity  vector  at  10-30  min  intervals,  and  moduli  of  this 
vector  averaged  over  [illegible]  min.  It  is  clear  that  such  instruments  are  not 
suitable  for  studying  the  "fine"  structure  of  the  velocity  field  and  high-frequency 
fluctuations  in  the  spectrum  of  sea  turbulence.  However,  the  extremely  large  hori- 
zontal dimensions  of  the  seas  and  oceans  cause  the  appreciable  large-scale  character 
of  horizontal  processes.  Therefore,  in  determining  the  large-scale  characteristics 
of  the  field  of  horizontal  velocity  components  in  the  ocean,  a long  series  of  obser- 
vations (on  the  order  of  many  days)  is  required,  and  short-period  fluctuations, 
which  make  a negligible  contribution  to  the  total  energy  of  horizontal  turbulence, 
may  be  neglected  altogether.  Therefore,  the  selection  of  BPV-2  should  be  considered 
fully  adequate  for  the  study  of  large  fluctuations  of  the  velocity  field  in  the 
ocean,  which  contain  a significant  fraction  of  the  energy  of  ocean  turbulence. 

Major  disadvantages  of  many  velocity  recordings  by  BPV-2  instruments  arc  the 
appreciable  errors  introduced  into  the  instrument  readings  by  their  mobility  on 
buoy-type  anchor  devices  used  in  the  measurement  of  sea  currents.  However,  as 
evidenced  by  experience  with  such  buoy  devices,  appreciable  errors  arc  observed  only 
in  the  readings  of  instruments  placed  in  the  upper  sea  levels  (approximately  down 
to  100  m).  Therefore,  we  selected  for  our  calculations  the  data  of  an  instrument 
which  operated  at  a depth  of  200  m,  where  the  errors  introduced  by  the  motions  of 
the  buoy  on  the  wave  may  be  considered  absent. 

Data  of  hydrological  current  meters  for  studying  large-scale  horizontal  tur- 
bulence in  the  sea  were  used  for  the  first  time  by  Shtokman  back  in  1939-40.  In 

1941,  Shtokman^  calculated  the  Reynolds  stresses  -pu'v',  -pu'2  and  -pv'2  for  sea 

turbulence.  This  was  done  by  using  data  of  day-to-day  current  measurements  with 

the  Ekman-Mertz  current  meter  in  the  Caspian  Sea.  The  measurements  were  made 

every  2 hours  and  lasted  23  days.  The  fluctuation  velocities  u'  and  v'  wore  deter-  / 1 7 

mined  by  Shtokman  by  subtracting  the  velocities  u and  v,  averaged  over  the  entire 

observation  period,  from  the  "instantaneous"  velocity  values  recorded  by  the  current 

2 3 

meter.  Later,  similar  studies  in  the  Black  Sea  were  made  by  Khlopov  and  Cezentsvey. 
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The  latter.  In  calculating  u'2  and  v'2,  successively  used  as  the  "Instantaneous" 
velocity  values  those  recorded  by  the  current  meter  and  those  averaged  over  periods 
of  1 h and  6 h.  The  fluctuation  velocities  were  obtained  by  subtracting  the  values 

of  the  velocity  u and  v,  averaged  over  the  entire  observation  period  (14  days),  from 
these  "Instantaneous"  values.  This  method  enabled  Gezentsvey  to  reveal  the  depen- 
dence of  u'2  and  v'2  on  the  averaging  period  of  "instantaneous"  velocities.  However, 
the  constancy  of  the  basic  averaging  scale  Tq,  used  in  the  calculation  of  u and  v, 

did  not  permit  her  to  relate  u'^  and  v'^  to  this  scale.  As  was  noted  above,  it  is 
precisely  this  dependence  which  should  be  of  Interest  to  an  Investigator. 

The  Reynolds  stresses  for  sea  turbulence  for  a constant  averaging  scale  Tn  were 

4 0 5 

also  calculated  by  Stommel  from  velocity  measurements  in  the  Straits  and  by  Ichlye, 

who  used  for  this  purpose  photographs  of  the  velocity  field  made  with  the  aid  of  an 

electromagnetic  current  meter  (EMCM)  in  the  region  of  the  Kuroshlo. 

The  quick-response  equipment  built  In  the  last  few  years  has  made  it  possible 
to  determine  Reynolds  stresses  caused  by  high-frequency  fluctuations  of  the  velocity 

field  in  the  ocean.  Studies  have  been  made  along  these  lines  by  Bowden  and  Fatrbalrn,^’ 

8—10  11  12  13 

Kolesnikov  et  al. , Grant,  Stewart  and  Moilliet,  ’ and  Bowden  and  Howe.  In 

these  authors'  experiments,  the  duration  of  the  current  velocity  recordings  usually 

did  not  exceed  5 to  10  min,  and  the  fluctuation  velocities  were  calculated  relative 

to  the  velocity  components  averaged  over  such  a time  interval.  The  majority  of 

these  studies  also  showed  that  for  such  scales  of  the  phenomenon,  the  statistical 

characteristics  of  sea  turbulence  obey  the  relationships  established  in  the  theory 

of  locally  isotropic  turbulence  by  A.  N.  Kolmogorov. 

In  studying  the  dependence  of  turbulent  Reynolds  stresses  in  the  sea  on  the 
averaging  scale  Tq,  it  is  natural  to  proceed  as  follows.  Let  the  recording  of  any 

given  component  of  the  current  velocity  at  some  point  of  the  sea  be  described  by  the 
function  u(t),  which  may  be  considered  continuous  and  differentiable.  Such  a func- 
tion may  be  represented  by  a superposition  of  simple  harmonic  oscillations  with 
different  amplitudes  and  periods  Tn  (or  angular  frequencies  U)q  = 211/T^) . If  by  virtue 

of  certain  consideration  we  choose  the  time  scale  of  the  averaging  equal  to  Tq,  we 
thereby  refer  all  components  with  periods  greater  than  Tq  to  "regular"  changes  in 

averaged  velocity  u(t),  and  higher-frequency  velocity  changes  should  in  this  case 
be  treated  as  turbulent  fluctuations  u'(t). 

Tn  order  to  separate,  in  a given  curve,  the  components  with  periods  greater 
than  Tq  from  other  oscillations,  it  is  obviously  necessary  to  "pass"  the  function 
u(t)  through  a certain  type  of  filter.  In  the  ideal  case,  low-frequency  oscilla- 
tions (with  <*>  < u)q  **  2tt/Tq)  should  pass  through  such  a filter  without  any  attenuation 

(the  amplitude  ratio  of  the  transmitted  to  the  arriving  signal  A/Aq  should  be  equal 

to  unity),  and  the  amplitudes  of  harmonics  witli  frequencies  greater  than  u>q  should 

fall  to  zero.  The  frequency  characteristic  of  such  an  "ideal"  filter  is  shown  in  / 1 71 
Fig.  1 by  a broken  line.  The  function  u(t)  can  be  filtered  either  by  means  of 
electronic  devices  (the  function  in  this  case  is  given  by  a suitable  electric  signal), 
or  by  analytical  means.  Analytically,  the  filtering  (or  smoothing)  operation  may  be 
represented  as  follows: 
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(2) 


5(t)  *(,)*, 

where  T is  the  integration  variable,  and  the  function  K(t)  is  called  the  smoothing 
core. 
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Fig.  1.  Spectral  characteristics  of  smoothing  cores  (3)  (broken  line)  and  (4) 
(solid  line). 


For  an  "ideal"  smoothing  with  a rectangular  frequency  (spectral)  characteristic, 
shown  in  Fig.  1 by  a broken  line,  the  function  K(t)  is 


K{x)  = 


sm  (»oT 
ax 


(3) 


where  » 2tt/Tq. 

Indeed,  if  such  a core  is  made  to  act  on  a simple  harmonic  oscillation  cos  uit, 
it  is  readily  seen  that  this  oscillation  remains  unchanged  when  o)  < Mg  and  will  be 
completely  "extinguished"  (the  integral  in  Eq.  (2)  will  be  zero)  when  u>  > (0Q.  When 
id  - nig,  the  integral  In  Eq.  (2)  will  take  on  the  value  of  1/2. 

In  a practical  application  of  the  smoothing  operation  to  experimental  data, 
the  use  of  an  ideal  filter  is  still  difficult  in  the  majority  of  cases.  Indeed, 
function  (3)  damps  out  relatively  weakly  with  increasing  T,  and  therefore  the  inte- 
gration in  formula  (2)  must  be  performed  over  a very  large  interval  of  change  in  T. 
Therefore,  instead  of  core  (3),  it  is  necessary  to  use  other  functions  different 
from  zero  only  on  a finite  segment  of  the  T axis.  In  this  case,  the  spectral  smoothing 
characteristic  deteriorates  to  some  extent,  and  becomes  different  from  the  "ideal"  one. 
On  the  advice  of  Professpr  W.  H.  Munk  (USA),  we  used  a smoothing  core  of  the  fol-  /I 
lowing  form,  very  convenient  for  practical  applications: 


K(x) 


1 -j-  COS 

- jr. 
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k _ _ n 

- — ^ T ^ | 

ID*  <i)# 


(4) 
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The  spectral  characteristic  of  the  smoothing  core  4 Is  shown  by  a solid  line 
in  Fig.  1.  It  Is  obvious  that  this  characteristic  differs  appreciably  from  the 
"ideal"  one.  When  core  4 is  used,  components  with  periods  greater  than  Tq  are  also 
partly  attenuated,  and  higher-frequency  components  do  not  damp  out  completely.  How- 
ever, as  is  evident  from  the  graph,  oscillations  with  frequencies  greater  than  2u>q 
will  be  practically  absent  in  the  smoothed  function,  since  the  characteristic  curve 
for  0)  > 2u)q,  oscillating  about  the  axis  of  abscissas,  approaches  it  rapidly  with 
increasing  oj,  and  even  the  first  maximum  has  an  absolute  value  of  less  than  0.03. 

At  the  point  o>  * U)q,  the  spectral  characteristic  of  core  4,  like  that  of  an  "ideal" 
core,  has  a value  of  1/2. 

As  already  noted,  the  initial  data  used  in  our  calculations  were  the  velocity 
recordings  made  by  the  BPV-2  instrument  on  a day-to-day  buoy  station  operated  in 
June-August  1958  in  the  Atlantic  Ocean  at  a point  with  coordinates  55*15'  N and 
16*30'  W and  a depth  of  3080  m.  A current  recorder  registered  the  velocity  vector 
modulus  and  its  direction  every  30  min.  These  data  were  used  to  calculate  the 
horizontal  velocity  components  u(t)  and  v(t),  which  were  then  subjected  to  a 
smoothing  operation.  The  smoothing  was  carried  out  with  core  4 using  values  of 
the  parameter  Tq  equal  to  3,  6,  12  and  24  h.  The  smoothing  technique  is  very 
simple.  First,  for  a chosen  parameter  Tq,  the  values  of  function  (4)  are  determined 
for  each  instant  XR  of  velocity  recording;  these  values  are  then  multiplied  by  the 
corresponding  values  of  u(tg  - x ) , and  the  products  are  summed.  The  sum  thus  cal- 
culated gives  the  value  of  u(tg)  for  a specified  value  of  time  tg.  Then,  the 
next  instant  tg  + At , where  At  - 30  min,  is  taken  as  the  initial  point,  and  the 
value  of  u(tg  + At)  is  similarly  obtained.  Continuing  this  operation,  we  eventually 

obtain  a series  of  discrete  values  of  the  function  u(t) , averaged  with  a given 
smoothing  period  Tq.  The  values  thus  obtained  for  smoothed  functions  of  horizontal 
velocity  components  are  given  in  Figs.  2 and  3.  In  each  case,  the  calculation  was 
made  for  200  consecutive  instants  tn,  i.e.,  the  interval  of  change  in  argument  t 
in  these  functions  is  100  h. 

The  graphs  of  Figs.  2 and  3 clearly  show  that  increasing  the  averaging  period 
causes  fluctuations  of  increasingly  larger  scale  to  disappear  from  the  curves,  and 
the  curves  themselves  become  increasingly  "smoother." 

Having  found  the  values  of  averaged  functions  u(t)  and  v(t),  one  can  also  easily 
obtain  the  values  of  the  corresponding  fluctuation  velocities  u'(t)  and  v'(t).  It 
is  sufficient  for  this  purpose  to  calculate  the  values  of  these  averaged  functions 
from  the  corresponding  "instantaneous"  values  of  u(t)  and  v(t).  Squaring  u'(t) 
and  v'(t)  (or  multiplying  them  out)  and  averaging,  we  obtain  the  values  of  Reynolds 
stresses  (we  neglect  the  minus  sign,  and  p = 1)  for  a given  averaging  scale  Tq. 

The  values  of  u'^,  v'^  and  u'v'  calculated  in  this  manner  for  Tq  « 3,  6,  12  and  24  h 
arc  given  in  the  table.  It  is  obvious  from  the  data  of  this  table  that,  as  expected, 
a regular  increase  in  the  turbulent  stresses  u'2  and  v1^  is  observed  with  increasing  /ly 

averaging  period  Tq.  At  the  same  time,  u'v'  remains  very  small  for  all  Tq.  This 

fact,  along  with  the  approximate  equality  of  u'^  and  v'^  (when  Tq  - 3,  6 and  12  h) , 
attests  to  the  isotropy  of  turbulent  fluctuations  of  these  dimensions  in  the  inves- 
tigated turbulent  current  in  the  ocean.  This  is  also  indicated  by  the  tabulated 
data  on  the  correlation  factor  between  u'  and  v' . Indeed,  this  factor  may  be  /171 

assumed  equal  to  zero  within  the  accuracy  of  the  calculation.  The  isotropy  of 
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Fig.  3.  Hie  function  v(t)  for  different  averaging  periods: 
a - "instantaneous"  values  of  the  function  v(t); 
b - with  averaging  periods  of  3 h; 
c - 6 h;  d - 12  h;  e - 24  h. 
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turbulent  fluctuations  breaks  down  to  some  extent  only  in  the  region  of  the  24-h 

ayeraglng  scale,  where  the  values  of  u'^  and  v'2  differ  appreciably.  This  phenomenon 
is  apparently  due  to  a strongly  manifested  singularity  in  the  region  of  the  velocity 
fluctuation  spectrum,  i.e.,  ordered  tidal  oscillations  of  the  diurnal  period,  which 
are  distinctly  different  in  the  graphs  of  Figs.  2 and  3. 

Data  on  the  magnitude  of  u'2  and  v'^  also  make  it  possible  to  analyze  the 
behavior  of  the  sum  of  these  quantities  with  increasing  Tq-  As  can  be  readily  seen, 
this  sum  is  merely  the  total  energy  E (per  1/2  unit  mass)  of  horizontal  turbulent 
fluctuations  with  periods  shorter  than  Tq.  Increasing  the  averaging  period  from  Tq 

to  T”  should  give  an  addition  to  the  turbulent  energy,  due  to  velocity  fluctuations 

with  periods  located  in  the  interval  AT  = T^  - Tq.*  If  the  smoothing  is  performed 

with  a whole  set  of  smoothing  parameters  Tg,  one  can  then,  generally  speaking,  plot 
a graph  of  the  function  E(T)  (or  E(w)).  In  other  words,  this  method  can  be  used  to 
study  the  energy  spectrum  of  large-scale  turbulent  fluctuations  of  the  velocity  field 
in  the  ocean.  If  however  the  increment  of  K (to)  due  to  turbulent  fluctuations  in 
the  frequency  interval  Aw  is  divided  by  Ato,  we  obtain  still  another  Important  charac- 
teristic of  turbulence  - the  spectral  energy  density  f(w)  = AE(<o)/Au).  Obviously,  in 


It  should  be  emphasized  that  the  statements  made  here  will  be  strictly  valid 
provided  an  "ideal"  smoothing  of  the  functions  u(t)  and  v(t)  is  used.  In  our  case, 
they  hold  only  with  a certain  approximation.  In  addition,  it  should  be  kept  in  mind 
that  a certain  error  is  also  introduced  into  the  values  of  the  function  E as  a result 
of  the  limited  duration  of  the  velocity  component  realizations  used  in  the  calculations. 


Fig.  4.  Experimental  values  of  the  Fig.  5.  Experimental  values  of  the 

function  E(w)  and  theoretical  curve  function  f(oj)  and  theoretical  curve 

of  the  2/3  power  law.  of  the  5/3  power  law. 

our  case,  only  three  values  of  the  function  f(w)  can  be  obtained,  which  were  found  / 1 7 1 
4 4 4 2 -1  * 

to  be  4.2-10  , 17. 3*10  , and  93.4‘10  cm  sec  . Since  the  A u)  values  used  for 
calculating  f (oo)  are  very  large,  it  remains  unclear  to  which  U)  one  should  refer  the 
values  of  the  function  f(io)  that  were  obtained  in  this  manner.  For  this  reason, 
the  graph  of  Fig.  5 does  not  show  points  of  the  function  f (to) , but  rather,  rectangles 
with  bases  equal  to  Aio,  and  heights  corresponding  to  the  values  of  f(w)  for  these 
intervals  of  change  in  the  variable  oj. 

The  values  of  the  function  E(u>)  are  shown  in  Fig.  4.  As  is  evident  from  Figs.  4 
and  5,  the  experimental  values  of  the  functions  E (oo)  and  f(oi)  are  in  good  agreement 
with  certain  regular  curves  whose  shape  will  be  specified  below. 

The  theory  of  locally  isotropic  turbulence  offers  expressions  for  the  functions 
E and  f which  hold  in  the  so-called  inertial  subinterval  of  turbulence  scales. 

These  expressions  are 

* * 

It  should  be  noted  that  the  functions  E and  f can  obviously  also  be  obtained  in 

the  usual  manner,  i.e.,  by  a Fourier  transformation  of  the  correlation  functions  for 
the  velocity  components  in  the  flow  in  question.  However,  our  method  of  smoothing  the 
functions  u(t)  and  v(t)  with  different  averaging  periods  Tq  has  the  advantage  that  it 
makes  it  possible  to  avoid  nonstationary  long-period  fluctuations  in  the  velocity  field 
and  in  addition,  along  with  the  values  of  the  functions  L and  [symbol  illegible],  one 

can  also  obtain  data  on  the  derivatives  3u/3t  and  3v/3t  for  different  Tq,  so  that  one 
can  in  turn  compare  the  orders  of  the  terms  in  the  Reynolds  equations,  as  will  be  dis- 
cussed in  detail  below. 
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E[k)  = ce’'-*-*/*,  f(k)  = cae’/'fc-*/-, 


(5) 


where  k Is  the  wavenumber;  E is  the  dissipation  rate  of  turbulent  energy  per  unit 
mass,  and  and  c ^ are  dimensionless  universal  constants. 

The  observations  whose  data  are  used  in  this  study  were  performed  at  a fixed 
point  of  space,  and  therefore  in  order  to  be  able  to  apply  formulas  (5)  to  these  / 1 7 1 

data,  it  is  necessary  somehow  to  replace  the  space  variables  in  Eq.  (5)  by  time 
variables.  This  substitution  may  be  made  with  the  aid  of  the  hypothesis  of  "frozen 
turbulence."  According  to  this  hypothesis,  turbulent  eddies  with  a space  wavenumber 
k during  their  motion  through  the  observation  point  give  rise  at  this  point  to  velo- 
city fluctuations  with  an  angular  frequency  0)  related  to  k by  k = 0)/V,  where  V is 
the  velocity  of  "transfer"  of  turbulent  eddies  by  the  mean  flow.  The  "frozen 
turbulence"  hypothesis  is  accurate  in  the  case  of  turbulent  flows  in  which  the  energy 
of  turbulent  fluctuations  is  much  lower  than  the  energy  of  mean  motion.  In  our  case, 
this  condition  is  not  fulfilled  very  well.  Indeed,  as  can  be  easily  seen  from  the 
tabulated  data,  the  rms  values  of  the  flow  velocities  at  high  values  of  Tq  approach 
the  mean  flow  velocities  (in  absolute  value).  This  fact,  along  with  the 
marked  variability  of  the  mean  velocity  with  time,  makes  it  necessary  to  apply  the 
"frozen  turbulence"  hypothesis  only  very  carefully  to  the  data  under  consideration. 
However,  we  will  use  this  hypothesis  anyway,  if  only  to  get  an  idea  of  the  possible 
form  of  the  dependence  of  E and  f on  the  frequency  u).  Formulas  (5)  are  thus  trans- 
formed to 

/?(o>)  = /(<o)  — a2&>~,/,1  (6) 

where 

at  — Cie’/,K’/>.  a2  — c2e,,lV’,t. 

The  graphs  of  functions  (6)  are  shown  in  Figs.  4 and  5.  The  values  of  coeffi- 
cients a^  and  a ^ were  chosen  as  follows:  a^  = 0.161  cm^  sec~8/3>  a ^ = 0.153  cm^  sec-®/^. 

As  is  evident  from  these  graphs,  the  experimental  points  are  in  very  good  agreement 

with  the  theoretical  expressions  (6),  with  the  exception  of  the  values  for  T_=  24  h 

-4  -1  u 

(oj  - 0.73" 10  sec  ),  which  are  located  slightly  above  the  theoretical  curves.  This 

may  again  be  explained  by  the  presence,  in  this  frequency  region,  of  tidal  currents 
giving  a definite  "spike"  in  the  spectral  functions  E(<o)  and  f(u)).  Obviously,  one 
might  suspect  that  such  good  quantitative  agreement  between  experiment  and  theory  is 
due  to  an  accidental  favorable  combination  of  errors  introduced  by  the  perturbations 
and  inaccuracies  listed  above.  However,  one  can  apparently  state  with  a fair  amount 
of  confidence  that  the  experimental  data  confirm  the  applicability  of  the  expressions 
of  the  theory  of  locally  Isotropic  turbulence  to  the  investigated  scale  interval  of 
horizontal  turbulent  formations  in  the  sea.  This  conclusion  should  not  be  particularly 
surprising  if  one  recalls  that  the  horizontal  dimensions  of  the  ocean  are  of  the  order 
of  thousands  of  kilometers,  whereas  the  spatial  dimensions  of  the  investigated  pheno- 
mena (calculated  by  multiplying  the  Tq  values  by  the  corresponding  moan  current  velo- 
city) are  only  1-1.0  km.  For  phenomena  with  such  spatial  dimensions,  the  hypothesis 
of  local  isotropy  had  also  been  confirmed  previously  by  checking  the  applicability 
of  the  4/3  power  law  to  the  coefficient  of  horizontal  turbulent  diffusion  in  the 
ocean. ^ 
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The  above  values  of  the  coefficients  a^  and  a ^ In  formula  (6)  make  It  possible 

to  estimate  the  order  of  magnitude  of  the  turbulent  energy  dissipation  rate  in  the 
sea,  t.  Considering  that  the  velocity  V is  of  the  order  of  10  cm/sec,  and  the  con- 
stants c and  c0  are  taken  to  be  equal  to  unity,  the  order  of  magnitude  of  e is 
1 *2-3 

found  to  be  10“^  cm  sec  , which  is  in  good  agreement  with  previous  estimates  of 
e in  the  sea.^^ 

The  last  two  rows  of  the  table  give  the  rms  values  of  the  time  derivatives  of  / 1 7 

of  the  velocity  components  u(t)  and  v(t)  averaged  with  a different  scale.  These 

derivatives  were  calculated  by  dividing  the  different  between  consecutive  values  of 

u(t)  and  v(t)  by  the  time  interval  separating  these  values.  It  follows  from  the 

table  that  the  rms  values  of  9u/3t  and  3v/3t  are  substantially  dependent  on  the 
averaging  scale,  decreasing  with  Tq.  This  result  illustrates  once  again  how  care- 
fully one  must  make  the  estimate  of  the  order  of  the  terms  in  the  equations  of  motion 
of  seawaters.  In  making  this  estimate,  it  is  always  necessary  to  specify  exactly 
the  scale  of  the  process  to  be  described  by  a given  system  of  equations  of  motion. 

An  analytical  expression  for  the  dependence  of  the  rms  value  of  3u/3t  (or  3v/3t) 

on  the  averaging  scale  Tq  may  be  obtained  from  the  following  considerations.  Obviously, 

3u/3t  is  determined  mainly  by  the  smallest-scale  velocity  oscillations  existing  in 

the  function  u(t)  (i.e.,  oscillations  with  period  Tq) . However,  according  to  the 

theory  of  locally  isotropic  turbulence,  the  amplitude  of  these  fluctuations  Au  is 
proportional  to  Tq  to  the  power  1/3.  And  since  the  characteristic  time  At  for  this 

fluctuation  is  again  equal  to  Tq,wc  obtain  for  the  rms  value  of  the  derivative  3u/3t 


where  Y is  a proportionality  coefficient. 

Expression  (7)  may  be  obtained  from  the  relation  between  the  rms  value  of  the 

derivative  3u/3t  averaged  over  the  interval  T and  the  structural  function  of  the 

velocity  field  D (x).  Indeed,  we  have 
uu 

1/(1  f*  )LmT. 

If  however  we  now  use  for  the  function  D (T)  the  2/3  power  law  known  in  the 

uu 

theory  of  locally  isotropic  turbulence,  we  immediately  arrive  at  Eq.  (7),  obtained 
above. 

A graph  of  the  power  law  (7)  with  proportionality  coefficient  y = 0.59  cm  sec 
is  given  in  Fig.  6.  It  Is  evident  from  this  figure  that  the  experimental  points 
are  again  in  good  agreement  with  the  theoretical  curve,  once  more  confirming  the 
applicability  of  the  theory  of  locally  isotropic  turbulence  to  horizontal  pro- 
cesses of  a given  scale  in  the  ocean. 
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Fig.  6.  Experimental  values  of  the  derivatives  3u/3t  (circles)  and  3v/3t  (crosses) 
for  different  averaging  periods,  and  theoretical  curve  of  the  2/3  power  law. 


The  values  obtained  for  the  nonstationary  terms  in  the  Reynolds  equations  of  / 17 1 
motion  may  be  compared  with  the  values  of  the  components  of  the  Coriolis  force, 
which  plays  an  important  part  in  the  dynamics  of  sea  currents.  Calculation  of  the 

_3 

components  of  Coriolis  accelerations  in  our  case  leads  to  values  ranging  from  1.8*10 

» -3  —2 

to  0.6*10  cm  sec  (depending  on  the  mean  velocity).  Comparison  of  these  figures 
with  values  of  nonstationary  terms  in  the  Reynolds  equations  leads  to  the  conclusion 
that  these  terms  surpass  the  components  of  the  Coriolis  force  only  for  averaging 
scales  (phenomena)  shorter  than  2-8  h.  For  larger-scale  processes,  terms  with 
the  Coriolis  force  become  predominant  in  the  equations  of  motion,  and  such  processes 
may  in  a certain  sense  be  considered  quasi-stationary. 

Tn  conclusion,  it  should  be  noted  that  the  numerical  estimates  given  in  this 
study  for  large-scale  characteristics  of  the  field  of  horizontal  velocity  components 
in  the  ocean  are  obviously  typical  only  of  the  dynamic  conditions  existing  during 
the  observation  period.  Under  different  conditions,  however,  such  estimates  may 
of  course  change  markedly.  However,  it  may  be  stated  with  sufficient  reason  that 
the  general  character  of  the  dependence  of  the  investigated  characteristics  on  the 
averaging  scale  Tq  is  fairly  universal  and  determined  solely  by  the  general  charac- 
teristics of  turbulent  liquid  motion. 
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